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Abstract 

We examine the use of distributions in numerical treatments of Anderson localisation and supply evidence that treating 
exponential localisation on Bethe lattices recovers the overall picture known from hypercubic lattices in 3d. 
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The first question when studying locaUsation of particles 
is how to characterise localised states in distinction to ex- 
tended states. Concepts as the localisation length or return 
probability are intuitive quantities, but restricted to non- 
interacting particles, and quantities as the conductance are 
hard to calculate, at least for disordered interacting systems. 

Recalling Anderson's first ideas about localisation [1] a 
characterisation using distributions of the local density of 
states (LDOS) pi as the basic quantity is appealing since 
such a description can be adopted to interacting systems as 
well. This leads us to consider two different approaches to 
the localisation problem given by the Anderson model H = 



with box distribution PUi) 



9(7/2 -|6,|)/7. 

The AAT-scheme [2] formulates the localisation problem 
directly in terms of distributions. It provides a selfconsis- 
tency equation for the distribution of the local (retarted) 
Greensfunction Ga which can be solved via Monte-Carlo- 
sampling. By its construction on a Bethe-lattice, being a 
loop-free approximation to hypercubic lattices, AAT can be 
understood as a kind of mean field approach to the locali- 
sation problem. Clearly, neglecting all nontrivial loops pre- 
vents any treatment of weak localisation. In contrast, expo- 
nential localisation can be treated to a large extent, as shall be 
demonstrated here. To compare with localisation on hyper- 
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Fig. 1. Distribution of LDOS for strong disorder (band center uj — 0.0, 
7 — 2.0 ^ 0.66 X 7crit). Arrows indicate the most probable value (mpv), 
the typical density of states ptyp, and the average density of states Pavc- 
The inset displays ptyp nearly at the localisation transition (uj — 0.0, 
7 — 3.0) versus the number of iterations Nuer in the Monte-Carlo 
procedure, for samples sizes 5x10, 10 , 10 from bottom to top. 
Increasing the sample size allows to separate localised and extended states 
ai'bitrarily close to the localisation transition. 

cubic lattices a refinement of the kernel polynomial method 
(KPM) serves as a quasi-exact numerical procedure to con- 
struct the LDOS-distribution without calculating wavefunc- 
tions of the disordered system, hence reducing both compu- 
tational time and memory demands (for details see [3]). 

In Fig. 1 the distribution of the LDOS for strong disorder 
is shown. The appearance of heavy tails and a high peak at 
small values indicates the qualitative behaviour of the dis- 
tribution at the onset of localisation. Increasing the disor- 
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Fig. 2. Mobility edge trajectory for the Anderson model calculated with 
KPM for a 3d-lattice (diamonds), and with AAT (circles), representing 
the distribution through a sample with 5 X 10 elements. The crosses 
indicate points in the {u. 7)-plane, which correspond to delocalised states 
if sampled with 2.5 X 10 elements. The dotted (dashed) line shows the 
mobility edge trajectory resulting from TMT (simplified TMT-variant). 

der, thus pushing the system over the localisation transition, 
causes the distribution to be singular, being entirely concen- 
trated at Pi = 0. Evidently the arithmetic mean pavo = (pi) 
is inappropriate to detect this qualitative change. But other 
quantities, e.g. the most probable value or the geometric 
mean ("typical density of states" ptyp = exp(lnpi)), that 
put sufficient weight on small values of pi, will be critical 
at the transition. We have discussed elsewhere [4] how to 
use the typical density of states in combination with a thor- 
ough investigation of its limiting behaviour close to the real 
energy axis to extract the position of the localisation transi- 
tion from the distribution (which allowed for a calculation 
of mobility edges in an interacting system). 

Fig. 2 shows mobility edge trajectories calculated with 
AAT and KPM. Both trajectories show the same qualitative 
features. In contrast to KPM the AAT scheme does not suf- 
fer from limitations due to the finite geometric size of the 
system under consideration. Though the resolution which 
the distribution is sampled with in AAT is restricted by the 
finite size of the Monte-Carlo sample, this limitation can be 
easily overcome by increasing the sample size (cf. Fig. 1, 
2). In principle this allows for an extremely precise deter- 
mination of mobility edges within AAT. 

Although both KPM and AAT construct the full distribu- 
tion of the LDOS the localisation transition already man- 
ifests itself in a suitable moment thereof (e.g. the typical 
density of states). As a matter of course, this triggers the 
question whether the localisation problem can be treated by 
an effective theory using this moment as an (scalar and lo- 
cal) order parameter. 

Indeed, the typical medium theory (TMT) [5] which is en- 
tirely based on the typical density of states is believed to cap- 
ture the most basic effects of Anderson localisation, and has 
been recently applied to the Anderson-Hubbard model [6]. 
From a pragmatic point of view it is quite favorable to re- 
place AAT by the simpler TMT. However, the typical den- 



sity of states is no obvious choice for the order parameter. In 
contrast to what is sometimes claimed it provides no good 
approximation to the most probable value (cf. Fig. 1). So we 
tried to simplify even further and replaced the typical density 
of states with the minimal values of the LDOS calculated for 
the maximal value of the on-site potential e^, i.e. ei = ±7/2 
for the box distribution. This means to incorporate only the 
maximal scattering contribution, originating from the deep- 
est impurity. We can then calculate a mobility edge trajec- 
tory which, at a first glance, does equally well compared 
with the one resulting from TMT (cf. Fig. 2). Nevertheless, 
the two approaches fail to obtain the reentrant behaviour of 
the mobility edges. It seems that both TMT and our sim- 
plified TMT-variant do only capture localisation caused by 
deep impurities ("pushing states outside the band"), but not 
the subtle effects caused by quantum interference. Consis- 
tent with this interpretation the critical disorder predicted by 
these approaches is smaller than given by AAT. Therefore 
we come to the conclusion that the treatment of disordered 
systems within TMT might be inherently problematic. 

To summarise: AAT is a theory of high precision for ex- 
ponential localisation on the Bethe lattice which recovers all 
qualitative effects of exponential localisation on hypercubic 
lattices. The use of distributions in numerical treatments of 
Anderson localisation is further justified by a comparison 
to results obtained for 3d-hypercubic lattices. While local- 
isation effects originating from local correlations are fully 
contained in AAT they are partially lost in averaged descrip- 
tions as TMT. 
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